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It is well known [5], pQ that in models with time-homogeneous local volatility functions and 
constant interest and dividend rates, the European Put prices are transformed into European Call 
prices by the simultaneous exchanges of the interest and dividend rates and of the strike and spot 
price of the underlying. This paper investigates such a Call Put duality for perpetual American 
options. It turns out that the perpetual American Put price is equal to the perpetual American 
Call price in a model where, in addition to the previous exchanges between the spot price and 
the strike and between the interest and dividend rates, the local volatility function is modified. 
We prove that equality of the dual volatility functions only holds in the standard Black-Scholes 
model with constant volatility. Thanks to these duality results, we design a theoretical calibration 
procedure of the local volatility function from the perpetual Call and Put prices for a fixed spot 
price xq. The knowledge of the Put (resp. Call) prices for all strikes enables to recover the local 
volatility function on the interval (0,xo) (resp. (xq,+o6)). 

Keywords: Perpetual American options, Dupire 's formula, Call-Put Duality, Calibra- 
tion of volatility, Optimal stopping. 



Introduction 

In a model with local volatility function q(t,x), interest rate r and dividend rate S 

jdSf = s(t, Sf)SfdW t + (r - 5)S?dt, t > 
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h(T,y)=E[e~ rT (y 
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of the European Put option considered as a function of the maturity T > and the Strike 
y > solves Dupire's partial differential equation [3] : 



d 2 yy h(T, y) + (8- r)yd y h{T, y) - 8h(T, y), T, y > 



d T h(T,y) 

h(0,y) = (y-x)+, y>0 

One easily deduces that the function h(T — t, y) for (t, y) G [0, T] x ]R^_ satisfies the pricing 
partial differential equation for the Call option with strike x and maturity T in the model 



y,T 



dS, 
6y,T 



?(T - t, Sf T )Sf T dW t + {5- r)S^ dt, t e [0, T] 



(2) 



with local volatility function q(T — t,y), interest rate 5 and dividend rate r. Therefore 



h(T,y) 



E 



8T/ny,T 



(5? 



r 



and one deduces the following Call-Put duality relation 



which is also a consequence of [T] 

VT > 0, Vx, y > 0, E [e- rT (?/ - S^)+] 



E 



x 



Since it derives from Dupire's formula, this Call-Put duality equality is closely related to 
calibration issues. One remarks that in the particular case of a time-homogeneous volatility 
function (q(t,x) = cr(x)), then Sf T also evolves according to the same time-homogeneous 
volatility function. 

In this work, we are interested in deriving such a Call-Put duality relation in the case of 
American options and in investigating consequences in terms of calibration. In the Black- 
Scholes model with constant volatility q(t, x) = ex, when r denotes a bounded stopping-time 
of the natural filtration of the Brownian motion (Wt)t>o, one has 



E 



E 
E 



X 



where the second equality follows from Girsanov theorem. Taking the supremum over all 
stopping-times r smaller than T one deduces that the price of the American Put option 
with maturity T is equal to the price of the American Call option with the same maturity 
up to the simultaneous exchange between the underlying spot price and the strike and 
between the interest and dividend rates. Extensions of this result when the underlying 
evolves according to the exponential of a Levy process have been obtained in [7j. Let 
us also mention that another kind of duality has been investigated in [TB]. But, to our 
knowledge, no study has been devoted to the case of models with local volatility functions 
like (HJ). 

In the present paper, we consider the case of perpetual (T = +oo) American options in 
models with time-homogeneous local volatility functions s(t,x) = cx(x). In the first part, 
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we recover well-known properties of the perpetual American call and put pricing functions 
by extending an approach recently developed by Beibel and Lerche [2] in the Black-Scholes 
case. This makes the paper self-contained. 

In the second part, we introduce the framework used in the remaining of the paper. 
In the third part of the paper, we consider the exercise boundaries as functions of the strike 
variable and characterize them as the unique solutions of some non-autonomous ordinary 
differential equations. 

The fourth part is dedicated to our main result. We prove that the perpetual American Put 
prices are equal to the perpetual American Call prices in a model where, in addition to the 
exchanges between the spot price of the underlying and the strike and between the interest 
and dividend rates, the volatility function is modified. We also derive an expression of 
this modified volatility function. Notice that in the European case presented above, time- 
homogeneous volatility functions are not modified. 

The fifth part addresses calibration issues. It turns out that for a given initial value xq > 
of the underlying one recovers the restriction of the time-homogeneous volatility function 
cr(x) to (0, xq) (resp. [x ,+oo)) from the perpetual Put (resp. Call) prices for all strikes. 
In the last part, we show that at least when 5 < r, in the class of volatility functions 
analytic in a neighbourhood of the origin, the only ones invariant by our duality result are 
the constants. This means that the case of the standard Black-Scholes model presented 
above is very specific. 

Acknowledgements. We thank Damien Lamberton (Univ. Marne-la-vallee) and 
Mihail Zervos (King's College) for interesting discussions. We also thank Alexander Schied 
(TU Berlin) for pointing out the work of Beibel and Lerche [2j to us and Antonino Zanette 
(University of Udine) for providing us with the routine that calculates American option 
prices. 

1 Perpetual American put and call pricing 

We consider a constant interest spot-rate r that is assumed to be nonnegative and an asset 
St which pays a constant dividend rate 5 > and is driven by a homogeneous volatility 
function a : — > that satisfies the following hypothesis. 

Hypothesis (7i vo i): a is continuous on M^j_ and there are < a < a < +oo such that: 

Vx > 0, g_ < a(x) < a. 

In other words, S t is assumed to follow under the risk-neutral measure the SDE: 

dS t = S t {{r - 5)dt + a(S t )dW t ). (3) 

With the assumption made on a, we know that for any initial condition x 6 R+, there is 
a unique solution in the sense of probability law (see for example Theorem 5.15 in [TU] . 
using a log transformation) denoted by (Sf,t > 0). Moreover, Theorem 4.20 ensures that 
the strong Markov property holds for (Sf,t > 0). Under that model, we denote by 

P a (x,y) = sup E [e- rT (y - ^) + ] and C a {x,y) = sup E [e~ rT {S x T - y) + ] 
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respectively the prices of the American perpetual put and call options with strike y > and 
spot x. Here, T oo simply denotes the set of the stopping times with respect to the natural 

filtration of (Sf,t > 0). Since e~ rt S? = xexp (-8t + J* a(S*)dW u - \ J*cr 2 (S*)du) and 
f a 2 (S*)du > a 2 t — > +oo, it follows from the Dubins-Schwarz theorem that 



-rt ox 



e~ rt S? -> Oa.s. (4) 

As a consequence, e~ rt (y — Sf) + — > li r= ^y and e~ rt (Sf — y) + — > Oa.s. On {r = oo}, 
we thus set 

e- rT (y - S*) + = l {r=0} y and e' rT (S* - y) + = 0. (5) 
Let us consider the second-order ordinary differential equation 

^a 2 (x)x 2 f"{x) + (r - S)xf(x) - rf(x) = 0, x > 0. (6) 

According to Borodin and Salminen ([3], chap. 2) the functions 

n , \ E[e-"f], if a? < 1 , , , , f 1/E[e-^], if x < 1 

VX > °> ^ = { 1/E[«A if x > 1 and ^ = { E[e-W], if x > I ( 7 ) 

where for x,y > 0, = mi{t > 0, 5*f = y} ( inf = +oo ), are the unique solutions (up 
to a multiplicative constant) that are positive and respectively increasing and decreasing. 
The volatility function a being continuous, ([7j) ensures that these functions are C 2 on M*,. 



Remark 1.1. It is easy using the strong Markov property to get 
Vx,y>0, E[e~ rT S} = 



fr( x )/fr(y), tfx<y 

fl(x)/fi(y), ifx>y. 



Remark 1.2. Assuming r > 0, one has (see Borodin and Salminen J3J/j 

lim Jl( x ) = lim fl( x ) = 

X— Hj x— >+oo 

lim/|(x) = 0, lim ft(x) = +oo. 

x^O x— >+oo 



T/je function /j (resp. zs i/ras, up to a multiplicative constant, the unique solution 
to (0|) suc/i i/iai lim /(a;) = (Vesp. lim /(a;) = 0/ 

x— ++oo x— >0 

Remark 1.3. In i/je case 5 = 0, we /jat>e £/ie analytical solutions : 



h( x ) 



if(x) 



where ip(x) = x 



f + °° ( 1 




/, b exp 


/ ua 2 {u) dU _ 







dv, f-\(x) = x. 
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Indeed, since f(x) = x is solution of ^a 2 (x)x 2 f"(x) + rxf'(x) — rf(x) = 0, we search 
a general solution that can be written fix) = xf(x). This leads to ^cr 2 (x)xf"(x) + (r + 
a 2 (x))f(x) = and then f'(x) = § exp [- /* J^du] . Therefore, 3C U C 2 G R, f(x) = 



C 2 x + C x x f+°° (± exp [- /; ^fadu 



dv. 



Now, we are in position to show the existence of an optimal stopping time and give the 
call and put prices. Let us mention here that the problem of perpetual optimal stopping 
is treated in the paper of Dayanik and Karatzas [4], for a general payoff function and an 
underlying evolving according to a general one-dimensional time homogeneous diffusion 
process. Villeneuve [14] considers a model where the constant dividend rate 5 in ([3]) is 
replaced by a function 5 (St) and gives sufficient condition on the payoff function ensuring 
that a threshold strategy is optimal. Here, we give a direct proof that generalizes the 
approach developed by Beibel and Lerche [2] in the Black-Scholes case. 

Theorem 1.4. Assume r > 0. For any strike y > 0, there is a unique x*(y) G (0, y)n(Q, ^y] 
such that = inf{t > 0,5^ < x* a (y)} (convention inf0 = +oo) is an optimal stopping 
time for the put and: 

Vx<x:(y),P a (x,y) = (y-x) + , Vx > x* a (y), P a (x, y) = ^Z^M/^z) > (y - x)+. (8) 
In addition, we have f[(x*(y)) < and: 

fi(K(y)) 

Last, the smooth-fit principle holds: d x P a (x* a (y),y) = — 1. 

Remark 1.5. Ifr = 0, Vx, y > 0, P a (x, y) = y since for any stopping time t, (y — S*) + < y 
and equality holds for r = +oo by (GJ). 



Proof. Let us define: 



(y - z ) + 

Wz > 0, h(z) = — — — — and h* = sup h(z). 

Jl( z ) z>0 



Since the function h is continuous such that h(y) = and h(0+) = (Remark 11.21) . 
x*(y) = sup{z > 0, h(z) = h*} belongs to (0, y) and is such that h(x*(y)) = h*. Since the 
function h is C 2 on (0, y), we have h'(x*(y)) = and h"(x*(y)) < 0. These conditions give 
easily 

h(xl(y)) + (y- xl(y))f{[xl(y)) = and f{ (x* a (y)) > 0. 

Since f± is positive and x*(y) < y, we have f[(x*(y)) < and deduce ©• The second 
order condition and equation ([6]) then give x*(y)(r — 5)f[(x* a (y)) — rfi(x*(y)) < and so 

r V - b%l(y) > 0- 
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Now let us check the optimality of and consider r G T 0<oo . By Fatou's lemma and 
Doob's optional sampling theorem, we have 



E[e— (y-S T )+] < liminfE [e"" At (y - S? At ) + ] 
= liminfE [e— "f^MS?, 



< /»(<(y))]imiiifE [e— M /|(^ At )] < Kx'MWx) 

since e~ rt fi(Sf) = fi(x) + f Q e~ ru a(S^)S^f'^(S^)dW u is a nonnegative local martingale and 
therefore a supermartingale. If x > x^(y), we have using Remark 11.11 

E[ e — "(y - S£,) + ] = Ele-^SMCy - 5^ ) + ] = (y - <(y))E[ e -^)] = h(&(y)Wx) 

and t£ is optimal for x > x*(y). Since x*(y) = sup{z > 0, h(z) = h*}, we have (y — x) + = 
h(x)fi(x) < fi(x)h(Xp(y)) for x > x^(y), and finally deduces (JHj) for x > x* a (y). 

We consider now the complementary case x G (0, x*(y)), and set r G 7o i00 - Using the 
strong Markov property and the optimality result when the initial spot is x*(y), we get 

E[e—(y - < E[e-^W(2/ - )+]. 

On {t < r*. (y) }, we have 5f < z*(y), d e -' rt (y - 5f ) = e^SS? - ry)dt - e- rt a(S?)S?dW t 

<o 

and so E[e" rTAr ^) (y - S* At * )+] < lim infE[e~ rTAT ^) A *(y - At )+] <{y-x). □ 

Now, we state the similar result for the call prices. 
Theorem 1.6. Assume 5 > 0. For any strike y > 0, there is a unique T*(y) G (y, oo) D 

and; 



?y, +oo) suc/i i/iat r^ 7 = inf{t > 0, > Y*(y)} zs an optimal stopping time for the call 



T*(v) - y 

Wx > r;(y),C a (x,y) = (x-y) + , Vx < Tl(y),C a (x,y) = f'pt^ M*) > ^~v) + ■ ( 10 ) 
In addition, we have /j(T*(y)) > and: 

T *«(y) - y = j^tty ( n ) 

Lasi, iae smooth-fit principle holds: d x C a (T*(y),y) = 1. 

Remark 1.7. If 5 = 0, Vx,y > 0,Co.(x, y) = x. Indeed, the Call-Put parity K[e~ rt (Sf — 
y) + ] = x — ye~ rt + E[e~ r *(y — azves the convergence to x in both cases r > 

and r = u>aen i — > +oo. iVow;, thanks to the Fatou lemma, we have for r G T 0tOO : 
E [e- rr (^ - y)+] < liminfE [e~ r ™{S* At - y) + ] < liminfE [e' rrAt S^ At \ = x. 

t^+oo t— »+oo 
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Proof. The proof works as for the put, and we just hint the differences. We define 

( Z - y) + 

\/z > 0, h(z) = — — — — and h* = sup h(z). 

J]\ z ) z>0 

Let us admit for a while that lim h(z) = 0. Then, since h(y) = and h is continuous, 

z— >+oo 

h reaches its maximum in T*(y) = inf{z > 0, h(z) = h*}, and T*(y) € (y, oo). This 
gives ffTTl) . We then consider the case x < T*(y) and show that t£ = Ty*, y s is optimal. 
Note that in the special case r = 0, we have to use Proposition 11.81 which is stated below. 
Finally, we prove that t£ is optimal when x > T* (y) using that <5Y* (y) — ry > 0. 

Now, let us check that lim h(z) = 0. In the case r = 0, it is straightforward using the 

2— > + 0O 

explicit form given in Proposition 11.81 below that f^(x) > 1 2s_ x 1+2S ^" — 1 for x > 1, and 
we have then lim h(z) = 0. When r > 0, Ito's Formula gives 

z— >+oo 



de- r \S l t ) 1+a = e- rt (Sl) 1+a {(a + l)a(S})dW t + [a(r + (a + l)a 2 (^)/2) - (a + l)5]dt} . 



When a > 0, the drift term is bounded from above by a(r + (a + l)cr 2 /2) — (a + 1)<5 and 
we can find a > such that this bound is negative since a(r + (a + l)o r2 /2) — (a + 1)5 — > 

—<5 < 0. Then, for x > 1, we have E[e~ rT:,:Ai (S'^i A4 ) 1+a ] < 1 thanks to Doob's optional 
sampling theorem. The Fatou lemma gives thenE[e-"*(^i) 1+a ] < 1, and therefore we get 
fl(x) = 1/E[e~ rr -] > x 1+a . This shows lim h(z) =0. □ 



Proposition 1.8. In the case r = 0, the unique nonincreasing and increasing solution 
of (0|) starting from 1 in 1 are respectively: 



= 1, / T (x) 



where ip(x) 





" r 25 d ~ 


/ exp 




/o 


J l ua 2 {u) 



dv. 



Moreover, we have f]{x) = 1/P(r^ < +oo) for x > 1 and = ¥(rf < +oo) /or 

are (0,1]. 



Proof. When r = 0, the differential equation \o~ {y)y f"(y) — 5yf'(y) 



is easy to 



integrate: f'(y) = C 3 exp 



y 25 



-du 



ua' 2 (u) 

g(y) = c 4 + c 3 



and then 



ry 


' r 26 du 


/ exp 




Jo 


J 1 ua 2 {u) 



dv 



for C 3l Ci G R. For x > 1, fi{Sli At ) is a bounded martingale that converges almost surely 
to /|(x)l{ r i <+00 } thanks to (jl]) and lim/^(x) = 0. Therefore fi(x)F(r^ < +oo) = 1, and 
the proof is the same for x e (0, 1]. □ 



Call-put duality for Perpetual American Options 



8 



To conclude this section, we state a comparison result which will enables us to compare 
x*(y) and T*(y) with the exercise boundaries obtained in the Black-Scholes model with 
constant volatility q_ (resp. a) where q_ (resp. a) bounds the function a from below (resp. 
above). 

Proposition 1.9. Let us consider two volatility functions 0\ and o<i such that Vx > 
0, o~i(x) < a 2 (x) and that satisfy (TL vo \)- We also assume that either /[' or /(' (resp. 
either f'J^ or f" U2 ) are nonnegative functions and r > (resp. 5 > 0). Then, we have 

Vx, y > 0, P ai (ar, y) < P CT2 (z, y) (resp. C ai (x, y) < C a , 2 (x, y)) 

and we can compare the exercise boundaries: 

Wy > 0,<(y) > x* a2 (y) (resp. T* ai (y) < T* a2 (y)). 

Here and in the proof below, we add in the notation for each mathematical object the 
volatility function to which it refers. El Karoui and al. [B] and Hobson [5] prove that for 
a convex payoff function, the price of an American option with finite maturity is a convex 
function of the underlying spot price. They deduce monotonicity with respect to the local 
volatility function. Their results imply at the same time the convexity assumption made 
in the above proposition and its conclusion. In this paper, we prefer to give autonomous 
proofs of these results in our simple framework. And we will first use proposition 11.91 to 
compare with the Black-Scholes case where convexity is obvious. We can then deduce 
(Lemma l3.ip that /(' and f" are positive for any a satisfying (TC vo i). 

Proof. Let us consider for example the put case with f" ai > 0. Let x > z > 0. Ito's 
formula gives: 



de- rt f UtTl (Sr 2 ) 



e- rt fi ai (sr 2 )°2(sr 2 )sr 2 dw t + e 



-rt 



2(nx,a 2 



dt. 



The term between brackets is nonnegative since we have G\ < cr 2 and f^ ai is a convex 

.)] > fi,<n(x) using the 



-ru n t\T, 



function solving ([6]). Therefore we get E[e" 
optional sampling theorem, where u n = inf{t > 0, Sf ,<T2 > n} A n. Since v n — > +oo 

a.s. 

and /| )(Tl (S^'°^ x,<r 2 ) is bounded by fi, ai (z), Lebesgue's dominated convergence theorem 

gives then E[e~ rTz ] > f^ ai (x)/fi )(T1 (z) and so E[e~ TTz ' } > E[e~ rTz ] using Remark [Lll 
The same conclusion holds when /j_ )(Ta is convex by estimating Efe"^' f ltCT2 (S x i%)}. We 



then get the result: if P ai (x,y) > (x — y) + , P ai (x,y) = (y — x* (y))E[e m h^~>j < (y — 

x.cro x.an 

x* ai (y))E[e x °i {y) ) = E[(y - S x l% )+e x ^ (y) ] < P a2 (x,y). Now we just observe that 

x % 1 (») 

{x > 0, P ai {x, y) > (x — y) + } C {x > 0, P a2 (x, y) > (x — y) + } and thus x* Ay) = inf{x > 
0,P ai (x,y) > (x-y) + ) > inf{x > 0,P a2 (x,y) > (x - y)+} = x* (y). □ 
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2 Framework and notations 

We will present in this section the framework that we will consider in all the paper. To 
clarify the duality, we will use names that implicitly refer either to the primal (or "real") 
world, or to the dual world. This denomination has no mathematical meaning since, as we 
will see, there are no difference between them. On the contrary, from a financial point of 
view, natural variables such as the interest rate, the dividend rate have their true meaning 
in the primal world, while in the dual world they interchange their role. This is the reason 
why we also name the primal world "real" world. 

The primal ("real") world 

The primal world is the framework we just have described. The spot interest rate r is 
constant and nonnegative, and St is an asset which pays a constant dividend rate 5 > 
and is driven by a homogeneous volatility function a : IR+ — > that satisfies (Ti. vo i) under 
a risk-neutral measure. The prices of the perpetual American put and call are respectively 
denoted by P a (x,y) and C a (x,y), and their exercise boundary by x* a and T*. 

The dual world 

In the dual world, 5 plays the role of the interest rate and r of the dividend rate; x plays the 
role of the strike and y is the spot value of the share. Let rj : R* + — > R* + be an homogeneous 
volatility function that is also assumed to satisfy (Ti. vo i)- We consider then (S t ,t > 0) the 
solution of dS t = S t ((5 — r)dt + rf(S t )dW t ) that starts from y at time 0. Under that model, 
we denote respectively by 

p v (y,x) = sup E [e~ 5T {x -^) + ] and c v (y,x) = sup E [e - ^^ - x) + ] 

the prices of the perpetual put and call with strike x > and spot y. We can define, as 
in the primal world, and g^ as the unique decreasing (non-increasing when 5 = 0) and 
increasing positive solution to 

]-rj 2 (x)x 2 g"(x) + (5- r)xg\x) - Sg(x) = 0, (12) 

and we name the exercise boundaries ^(x) < x and y^(x) > x that are respectively 
associated with the put p v {y,x) and the call c v {y,x). 

Notations 

The aim of this paper is to put in evidence a duality relation and interpret put (resp. 
call) prices in the primal world as call (resp. put) prices in the dual world for a specific 
volatility function r\ = a (resp. r} — a). When r = (resp. 5 = 0), this is trivial since 
P a (x,y) = c v (y,x) = y (resp. C a (x,y) = p v (y,x) = x) but not really fruitful, and we take 
thus the following convention in the sequel. 
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Convention 2.1. We will always assume r > (resp. 5 > 0) to state properties on P a 
and c v (resp. C a andp n ). 

Both worlds being mathematically equivalent, we will work with the put price in the 
primal world and the call price in the dual world in order not to do the things twice. 
Following the Convention 12. 11 we will consider a positive interest rate. Let us then denote 
from now: 



/ = fi and g = g h 



and define a(y) = V f^*{y)) an< ^ @( x ^ = g(y* (x)) • ^ ne functions a and (3 are positive functions 
and it follows from the present section that : 

Vy>0, Vx>x*M, P a (x,y) = a(y)f(x) (13) 
Vx>0,Vy<y;(x), cr,(y,x) = P(x)g(y). (14) 

That product form will play an important role for the duality. Let us finally introduce 
notations relative the the Black-Scholes model. We define for q > 

, , 5 - r + q 2 /2 - s /{5-r + q 2 /2) 2 + 2r^ 2 n 
a\S) = ~ 2 < 0, 

, r-6 + e72 + y /(5-r- q 2 /2) 2 + 25q 2 , . . 1 

and we can easily check that f(x) = x a ^ (resp. g = x b ^) when u{x) = q (resp. r)(x) = q). 
In that case, the unique solution to @ (resp. ffTTj) ) is: 

= a(q) - l V ^ reSP ' = bU) - 1 ^' ^ 



With Proposition II. 9[ we deduce very useful estimations on the exercise boundaries. 
Lemma 2.2. If, \/x > 0, a < a(x) < a and q_ < r)(x) < a, then we have : 

-y < x* a (y) < ;~ 7 y with \~ J < mm(l,r/(5) (16) 



a(cx) — 1 a a (a) — 1 a(g_) — 1 

^ -x < yl(x) < -j^ - x uritfi max(l,g/r) < ,. 6 ^ - . (17) 



" ' 6(2) 

Proof. It is straightforward from Proposition II .91 and ( fl"5l) to get the bound on the exercise 
boundaries since x i— > x 1 is convex for 7 G]0, 1[. We have to show that for q > 0, J^)-i < 
min(l, r/8) and > niin(l, 8/r). Since a(c) < and b(q) > 1, we get that a(q)/(a(q) — 

1) G (0, 1) and b(q)/(b(q) — 1) G (1, +00). We can also check that a(<j)/(a(<?) — 1) is a root of 
the polynomial Q(X) = 5X 2 -(r+5+q 2 /2)X+r. AsQ(x)=0 a 2 x = 2(l-x)(r-5x) 

and since a(q)/(a(q) — 1) G (0, 1), we then deduce that a(<?)/(a(<?) — 1) < r/8. In the same 
way, we have b(q)/(b(q) — 1) > 8/r. □ 
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3 ODE for the exercise boundary 

We have seen previously that the exercise boundaries satisfy 

x*M - y = f(**M)lf(x*M) ( res P- y* v ( x ) ~ x = 9(yfa))/g'(y*(x))). (is) 

We will soon prove that for fixed y > (resp. x > 0), ( {TBI admits a unique solution 
x*(y) (resp. 

Lemma 3.1. TTie function f (resp. g') is negative (resp. positive) and f" (resp. g" ) is 
positive on (0, +oo). Moreover, the boundaries x*(y) and y*(x) are respectively the unique 
solution to y — x + f(x)/f'(x) = and y — x — g(x)/g'(x) = 0. Last, x*(y), a(y), yt{x) 
and /3(x) are C 1 functions on . 

Remark 3.2. Positivity of f" and g" and ( flgj) and [L~$\ ) imply positivity of d^.P (T {x,y) and 
dyC v (y,x) in the continuation regions. 

Differentiating ( fl8l) with respect to y (resp. x), one obtains 1 = (x*)'(y) ^ x °p(x*(y))^ V ^ 
(resp. 1 = (y*)\y) 9{v ^y^ iv)) ). Using and equation (Ej) (resp. (ED) one deduces 
the following result (see equation (|2"2"|) below). 

Proposition 3.3. Lei us assume that the volatility functions a and g satisfy (7i V oi)- Then, 
the boundaries x*(y) and yt(x) satisfy the following ODEs: 

(x *y (v ) = **MMK(y)) 2 (19) 



v 2 (y*r,(x))y*r,(x) 2 

' 2(y*(x) — x)(ry*(x) — 5x) 



Proof of Lemma \3.1\ We only give the proof in the put case, the argument being similar 
for the call. By ([6]), for x > 0, f"(x) has the same sign as h(x) = rf(x) + (5 — r)xf'(x). 
If for some x > 0, f'(x) = 0, then since / is positive, f"(x) > 0. Therefore x is a local 
minimum point of / which contradicts the decreasing property of this function. Hence /' 
is a negative function. 

When 5 < r, h and therefore /" are positive functions. When 5 > r, we remark that if 
f"(x) = then h'(x) = Sf'(x) < 0. Since the continuous function /" and h have the same 
sign, this implies that 

Vx > inf{^ > : f"{z) < 0}, f"(x) < 0. (21) 
Now for y > 0, by (jSJ) then ( [181) , we have 

x* a ( y ) 2 a(K(y)) 2 f"(K(y)) f&M) 



f>(x*(y)) f'(x*M) 



{r-8)x*M=5x*M-ry. (22) 
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By (fTBj) . the right-hand-side is negative and moreover lim y ^ +00 x* (y) = +00. Hence 
sup{;z > : f"(z) > 0} = +00 and with §ZB), we conclude that /" > 0. 
According to (TTBI) . F(x*(y),y) = where 

F(x,y) = y-x + f{x)/f{x). 

The function F is C 1 on (0, +00) x (0, +00) and such that 

Wx,y > 0, d x F{x,y) = -f(x)f" (x)//'(x) 2 < 0. 

Therefore for fixed y > 0, x*(y) is the unique solution to F(x, y) = 0. Moreover, y — » x*(y) 
is C 1 by the implicit function theorem. Last, one deduces from (|T8l) that a(?/) is a C 1 
function. □ 

The positivity of /" and g" gives the following result. 

Corollary 3.4. The comparison result stated in Proposition \1.9\ holds for any G\ < o 2 
satisfying (H vo i) . 

Let us now give a uniqueness result for the ODEs (|T9|) and (pOl) . 

Proposition 3.5. There is only one solution x* of (T72J] (resp. y* of ffify ) defined on 
(0, +00) that satisfies > 0, C\y < x*(y) < c 2 y with < c\ < c 2 < min(l,r/5) (resp. 
Vx > 0, d\X < yl{x) < d 2 x with d\ > max(l, 8/r) ). 

Proof. Let us first remark that the uniqueness result for ( |T9i) is equivalent to the uniqueness 
result for (|2Q|) . Indeed, it is easy to see that x* a {y) is solution of (|T9|) if and only if 
y(x) := l/(x*(l/x)) is solution of (|2"U]) with the volatility function 77 (x) = cr(l/x). This 
new volatility also satisfies (7i V oi)- Moreover, d\X < y(x) < d 2 x with di > max(l,5/r) if, 
and only if < c v y < x*(y) < c 2 y with < c\ < c 2 < min(l, r/5). 

Let us suppose then that there are two solutions of (1201) . yi(x) and y 2 (x), that are 
defined on M + and satisfy d 2 x > yj(x) > dix for some d 2 > d x > max(l,5/r). Since 
y[(x) > for x > 0, yi is invertible and we have: 

— y' 1 (y 2 (x)) = y2{xfr](y 2 {x)) 2 2(y 2 (x) - y^ 1 {y 2 (x)))(ry 2 (x) - Sy^jy^x))) 

dx 1 2(y 2 (x) — x){ry 2 {x) — 5x) y 2 (x) 2 r](y 2 (x)) 2 

Maj ~ Vl' 1 {y2(x)))(ry 2 (x) - 6y^ 1 (y 2 {x)) 
(y 2 (x) - x)(ry 2 {x) - Sx) 

Thus, the function i/j(x) = y : [ 1 (y 2 ( y x))/x solves 



ip'ix) = — 

X 

1 

X 



Uiix) — ip(x)x ry 2 (x) — Sxipix) 

T~\ x T~\ x ^ x 

y 2 [x) — x ry 2 [x) — ox 

ip(x) — 1 \ / ip(x) — 1 



y 2 (x)/x-lj\ ry 2 (x)/(5x) 



ip(x) 



(23) 
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The estimation d 2 x > Dj{x) > d\X for j G {1, 2} with di > max(l, 5/r) implies that : 

3A > 0, Vx > 0, l/A < if}(x) < A, (24) 

Vx > 0, ipix) < mm , — , 1 > and — 1 > 0. (25) 

\ x ox J x ox 

Since local uniqueness holds for ( 123]) by the Cauchy Lipschitz theorem, the only solution 
<p such that ip{l) = 1 is the constant if = 1. Therefore checking that ( 1241 does not hold 
for solutions <p satisfying (1251) and such that <^(1) 7^ 1 is enough to conclude that ip = 1. 

Let 90 be a solution to (!2"3"j) satisfying (I2"5"]) . If v?(l) > 1, by local uniqueness for (!2"3l) . 
for all x G M% y?(x) > 1. By ([25]), one deduces that for all x G M;, y/(x) < ±=£& < 0. 
Therefore, f'(x) < (1 — <p(l))/x for x G (0, 1], and we have 

cp(x) > <p(l) + (1 - <p(l)) ln(x) -> +00 

which is contradictory to ( 1241) . In the same manner, if y?(l) < 1, tp(x) < 1 for x G 
and </? is strictly increasing. In particular, for x < 1, v 9 '^) — (1 ~~ l P(^))/ x an d therefore 
(^(1) — y?(x) > (1 — y?(l)) ln(l/x) — > +00 and this yields another contradiction. □ 

x— >0 

Corollary 3.6. Let us denote C = {f G s.t /(0) = 0, 30 < a < 6,Vx > 

0, a < /'(^) < b}. The application o 1— > x* (Vesp. 1— ► 7/*^) is one-to-one between the 
set {a G C(R+) £/iat satisfies (7i V oi)} « n <^ se£ of function C x = {x G C, s.i 30 < ci < 
c 2 < min(l,r/5), Vy > 0, cyy < x(y) < C2y} (resp. C y = {y G C, s.t. 3max(l,5/r) < d\ < 
d 2 ,\/x > 0, dix < y(x) < d 2 x}.) 

Proof. If a is a continuous function satisfying (7i V oi), by ( fT9l and (fl6l) . x* belongs to C x . 
The one to one property is easy to get. If x* = x* 2 with o\ and o 2 satisfying (7i V oi), the 
ODE ([19]) ensures that af(x* ai (yj) = cri(x* 2 (y)) for ?/ > 0. Therefore a 1 = a 2 . 
Let us check the onto property and consider x*(y) G C x . The function a defined by 

a f x *f y \\ = \/ 2 (y ~ x*(y))(ry - 5x*(y))x*'(y) 

x*(y) 

is well defined thanks to the hypothesis made on x*. As x* satisfies ( fTT)]) and solves the 
same ODE ( fl9]) as x*, we have x* = x* using Proposition 13.51 

The proof for 77 1 — ^ ?/* is the same and gives incidentally the expression of 77 in function of 
the exercise boundary y*(x): 

v / 2{y*(x) - x){ry*(x) - 5x)y*'{x) 

viv 0*0) = 7r~\ • 27 

□ 
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4 The call-put duality 

This section is devoted to the key result of the paper : for related local volatility functions 
cr and r), we can interpret a put price in the primal world as a call price in the dual world. 

4.1 The main result 

Theorem 4.1 (Duality). The following conditions are equivalent: 
1. 

Vx, y > 0, P a (x, y) = c v (y, x). (28) 

2. x* and y* are reciprocal functions: Vx > 0, x*(y*(x)) = x. 

3. 7] = a where 

~ ( x 2(y - xl(y))(ry - 5x* a (y)) 
yx* a {y)a{x* a (y)) 



4. a = rj where 



2(y*(x) -x)(ry*(x) - Sx) 

T](x) = — - — . ' ., . (30) 

y*{x)xT](y*(x)) 



Remark 4.2. Thanks to relation [Tb}) (resp. [P7\)), if a (resp. rj) satisfies (7i V oi) then the 
dual volatility function a defined by / TffPj) (resp. rj defined by l(3U\) ) satisfies (7Y vo i). 

Proof. 1 2 : We have on the one hand P a (x,y) = y — x on {(x, y), x < x*(y)} and 

P a (x,y) > y — x outside, and on the other hand c v (y,x) = y — x on {(x,y), y > y^(x)} 
and c^iy.x) > y — x outside. The duality relation (|2"H]) imposes then that {(x,y), x < 

K(y)} = {&y), y> vfa)} and so v&*M) = y- 

2 => 3, 4 : Taking the derivative of the last relation, we get thanks to ( fTQl) and ( !20l) 
<(y) 2 °«(y)f n 2 (y)y 2 _ 1 1 deduce (12011 and (ROD 

2(y-xUy))(ry-8x*{y))2(y-x*(y))(ry-8x*{y)) ~ 1 ^ aeClUCe aI1Cl «' 

3^2 (resp. 4 2 ) : By (fT9]) (resp. (T20D) and ([20]) (resp. (p])), a;* -1 (resp. i/*" 1 ) 
satisfies (T5U]) (resp. (USD). Since by ([11]) (resp. (fT7]) ) this function satisfies flT7]) (resp. (PToT)). 
one concludes by Proposition 13.51 

2 =>- 1 : The equality ( |28|) is clear in the exercise region since {(x,y),x < x*(y)} = 
{(x,y),y > y*(x)}. Let us check that it also holds in the continuation region. Using the 
product form (fl4l) . and the smooth- fit principle (Theorem 11.61) we get for all y e 

y - <{y) = P(K(y))g(y) 

1 = -P(x*(y))g'(y). 

Differentiating the first equality with respect to y, one gets 1— x* (y)' = x* {%))' f3\x* a (y))g{y)+ 
f3{x* a {y))g'{y), which combined with the second equality gives 



-l = {3'(x:(y))g(y). 
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Dividing by the first equality and using (ITS]) , one deduces ^(x*(y)) = y(x*(|/)). Since 
x* : IR^_ — > is a bijection, there is a constant C ^ such that (3 = C f . Since 
\/y > 0,a(y)f(x* a (y)) = y - x* a (y) = P(x*(y))g(y), one has a = g/C. From (HBD and (HD, 
one concludes that (1281) holds. □ 



In this proof, we have shown that a is proportional to g, and so there is a constant 
C > such that P a (x,y) = Cf(x)g(y) for x < x*(y). In the Black-Scholes' case, we have 
P a (x,y) = C 'x a ^ 'y b ^ for x < x*(y), and we are able to calculate C using the boundary 
condition P a (x*(y),y) =y- x* a (y). We get that C = b(a)- b{a) / \-a(a)) a{a) . We retrieve 
then the already known analytical formulae for the put prices (e.g. Gerber and Shiu [8]). 

4.2 An analytic example of dual volatility functions 

By flHUl) and (I2U|) . if y* G C y (where C y is defined in Corollary 13. 6p . then the reciprocal 
function of y* is the put exercise boundary x* associated to the local volatility function 



(XX 



y/2(ry*(x) — 5x)(y*(x) — x) 



x^Jy*(x)' 

Now by ( 1271) . y* is the call exercise boundary associated with the dual volatility function : 



= y/2(y - x*{y))(ry - gggjgggg) 

y 

Let us consider the family of exercise boundaries 

*/ \ x + a 
y i x ) = X T — — 
ox + c 

where a, b, c are positive constants such that max(c/a, b) < min(l, r/8) (condition ensuring 
y* G C y ). Since y*(x)' = (bx 2 + 2cx + ac)/(bx + c) 2 , one has 



, ((r — 5b)x + ra — 5c)((l — b)x + a — c) 
a[x) — \l 2 — ,x > U. 

te z + 2cx + ac 

Moreover, the function x* a (y) is the only positive root of the polynomial function: X 2 
X(a — by) — cy, that is: 



and 



x l(y) = ^{by ~ a + V(°y - a ) 2 + A °y) 



y/% - x*(y))(ry - Sx* a (y))(bx* a (y) 2 + 2cx* a (y) + ac) 
My > 0, a(y) = 



y(bx* a (y) + c) 

This example enables us to check numerically the duality. We have plotted in Figured], the 
prices of an American put P a (T,x,y) in the primal world for the local volatility a{x) and 
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Figure 1: P a (T, x, y) and c&(T, y, x), and C CT (T, x, y) and pa-(T, y, x) as functions of the maturity 
T for x = 0.5, y = 0.4, r = 0.2, 5 = 0.1 and the volatility parameters (a,b,c) = (1,0.4,0.1). 

an American call Ca{T,y,x) in the dual world for the local volatility a(x) as functions of 
the maturity T. These prices have been computed using a finite difference method. We can 
see at T = 10 that the limit value is quite reached and both prices are equal. The plots are 
nonetheless distinct which means that the duality does not hold for finite maturities. We 
have also plotted, in function of T, C a (T, x, y) in the primal world and pg-(T, y, x) in the dual 
world to check numerically whether the volatility function a is such as C a (x, y) = p&(y, x). 
As we can see, the curves do not seem to converge toward the same limit when T is large. 
This means that the volatility function a such that Vx, y > 0, C a (x, y) = p&(y, x) (obtained 
from a as rj is obtained from 7] but with exchange of r and 5) is different from a. 

5 Consequence of the duality : A (theoretical) method 
of calibration for the volatility a(x) 

In that section, we will put in evidence the importance of the duality within the calibration 
scope. We suppose for this that we are on a (virtual) market where are traded perpetual 
securities, and where the short interest rate r and the dividend rate 5 can be observed. 
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This means that we know the price of the share Xq, and the market quotes on that share 
the perpetual American puts and calls for all strikes K > 0. We name respectively p(K) 
and c(K) these prices and denote: 

X = sup{K > 0, c(K) = x -K} and Y = M{K > 0, p{K) = K - x }. (31) 

We will first suppose that the put and call prices derive from a time-homogeneous local 
volatility model before relaxing this assumption. 



5.1 The calibration procedure 

Let us assume that there is a volatility function a satisfying (H vo i) such that for all K > 0, 
p(K) = P a (xo, K) and c(K) = C a (xo, K). The following proposition says that these prices 
characterize a and its proof gives a constructive way to retrieve the volatility function from 
the prices. 

Proposition 5.1. Let us consider xq > 0. The map 

a h- ((P a (x , K), C a (x , K)),K> 0) 
is one-to-one on the set of volatility functions satisfying (7Y vo i)- 

Proof. We first consider the put case. The differential equation satisfied by the put prices 
in the continuation region makes only appear the values and the derivatives in x, K being 
fixed. Hence, we cannot exploit directly the prices. But the duality relation enables to 
get a differential equation in the strike variable. Thanks to the Duality Theorem, we 
have P a (x ,K) = c„(K,xq) for some a satisfying (7Y vo i)- It is then easy to calibrate er(.). 
Indeed, one has K **f ? p"{K) + K{5 - r)p'(K) - 8p(K) = for K < Y = yi(x ). Since 
the differential equation is valid only for K < Y , we only get a on (0, Y] by continuity: 



v^/v -fin 1 Wp(K) + K(r-5)pr(K)) 

which is well defined since p"{K) = d\ca{K,x) > (Remark 13.21) . Then, we can calculate 
the exercise boundary for x G (0,x ], solving (!20l) supplemented with the final 

condition y%{xo) = Y backward. This step only requires the knowledge of a only on the 
interval (0, Y]. Finally, we can recover the desired volatility o~(x) for x < Xq thanks to 



^^^■ ^^ar 1 (32) 

Now let us consider the calibration to the call prices. This relies on the same principle, 
but we have to be careful because the Duality Theorem is stated given to the call interest 
rate S and dividend rate r. So we have to interchange these variables when we apply 
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that theorem. There is a function d satisfying (Tt vo i) such that : WK > 0, C a (x ,K) 
Pz(K,x ). We have 

^K 2 d(x) 2 c"(K) + (5 - r)Kc'(K) - 8c{K) = 
for K > X = Q(x ). Thus, we get 

<2(5c(K) +K(r-6)d(K)) 




c"(K) 



which is well defined for analogous reasons. We can then obtain as before the exercise 
boundary solving ( fT9l) forward 



and we finally get the volatility a(y) for y > xq using the Duality Theorem. More precisely, 
we interchange r and 6 in ( l29i) to get 

yy) y&(y)m(y)) 

□ 

This calibration method, although being theoretical, sheds light on a striking and in- 
teresting result: the perpetual American put prices only give the restriction of a(x) to 
(0, x ] and the call prices only the restriction of <r(x) to [x , +oo). This has the following 
economical interpretation : long-term American put prices mainly give information on the 
downward volatility while long-term American call prices give information on the upward 
volatility. This dichotomy is remarkable. In comparison, according to Dupire's formula [5], 
there is no such phenomenon for European options : the knowledge of the call prices gives 
the whole local volatility surface, not only one part. In other words, the European call and 
put prices give the same information on the volatility while the perpetual American call 
and put prices give complementary information. 

Thus, one may think that the perpetual American call and put prices only depend on a 
part of the volatility curve. This is precised by the Proposition below that gives necessary 
and sufficient conditions on the volatility functions to observe the same put prices (resp. 
call prices). 

Proposition 5.2. Let us consider x > and cri(.) ; o- 2 {.) two volatility functions satisfying 
(H\o\)- Then, the following properties are equivalent: 

(i) \fy> 0, P ai (x ,y) = P a2 (x ,y) (resp. \/y > 0, C ai (x ,y) = C a2 (x ,y)) 

(ii) < y^ 2 (xo), <5"i(y) = 0^2 (y). (resp. Vx > £*"i( x o), = &2(x) where dj denotes the 
local volatility function such that Vx, y > 0, C a . (x, y) = p^ (y, x) .) 
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(Hi) Vx G (0,x ], (J\{x) = cr 2 (x) and z/z (xo) 
a 2 (x) andQ^xo) = £~ 2 (x ).) 

4>) 



(iv) \/x G (0,xq], a"i(x) = 02(2;) and 



cr 2 (x) and 



/t^Oo) /t,o- 2 (*o) 



2W * 



n| 2 (x ). fresp. Vie [x ,+oo), <7i(x) 
/I'^jao) - ^ res ^- Vx G [a?o,+oo), (7i(a;) 



^ f\ ai and fi tlJ2 (resp. f^ ai and f^ a2 ) are proportional on (0, x ] (resp. [x ,+oo)). 
(vi) Vx < x ,\/y > 0, P ci {x,y) = P a2 (x,y) (resp. Vx > x ,Vn > 0, C ai (x,y) = C a2 (x,y)). 

Remark 5.3. • Among these many conditions, let us remark that condition (ii) on the 
dual volatility is much simpler than condition (Hi) on the primal volatility since the 
latter requires the equality of the dual exercise boundaries at xq. 

• When 5 = 0, according to Remark \1.3l in the put case, condition (iv) also writes 
Vx G (0,x ], <J\{x) = cr 2 (x) and 



— exp 

v z 



2r 



uariu 



-du 



dv 



— exp 

v z 



2r 



uaz\u) 



-du 



dv. 



• Since, by definition of f^ >(r . (resp. f\ )0 ) and the strong Markov property V0 < z < 
x, E[e~ rT ^} = f Uaj (x)/fl aj (z) (resp.\o < x < z, E^V] = f\^(x) / f\^(z)), 
the probabilistic counterpart of assertion (v) is V0 < z < x < xo, E[e~' rr<T i' z ] = 
E[e- rT ^] (resp. Vx < x < z, E[e _rT *i-*] = E[e -r73 s.«]J. 

Proof. We consider for example the put case. 

(i) => (ii) : See the proof of Theorem 15. 11 

(ii) => (Hi) : Let us define ip(x) = (y^ 1 )~ 1 (y^ 2 (x))/x. We can show as in the proof of 
Proposition 13.51 that ip(xo) = 1 and then ip = 1 on (0, xq], otherwise it would go to or +00 
when x — > 0, which is not possible thanks to (IT7|) . We get then Vx G (0, x ), o- 1 (x) = a 2 (x) 
using flHUj) that express Oj in function of y~ . and dj, j G {1,2}. 

(Hi) => (iv) : Thanks to ( JTH1) and Theorem 14.11 we have °^ — ~ 1 

-1 = f'l,<r 2 ( X °) 

(iv) => (v) : 



/|,<7lO0) 



Vs x (aH))-aH) 



=> (v) : The set of solutions to -^a\(x)x f" (x) + (r — 5)xf'(x) — rf(x) = 
on (0, x Q ] is a two-dimensional vector space, but thanks to the relation jp 1 ^ — J^J^ji 
fl j(T1 and /| jCT2 are proportional on (0, x ] : 



Vx < X , fl,ax(x) = { i,fT1 ^°\ fi,a 2 (x). 



(V) (VI) 

then (y~ (x) — x) 



fl,(T 2 ( X 

The proportionality implies that Vx G (0, x ], 4^ 

Jl, (71 



O) fl,<T 2 {x) 

(y% (x) — x) _1 using (|18p and Theorem 14. 11 Therefore 



(34) 



and 



Vx G (0,x ], yUx) = y~Jx). 
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We have cn cri (y~ l (x))fi tai (x) = a a2 (y~ 2 (x)) f \,a 2 {x) using ([T3l . and obtain from ([Ml) that 

Vx < * , < vKxo), a ai (y) = f -j^\ acT2 {y) = fy^f\ a<T2 {y) . (35) 

Thus, we deduce from ( TTBl . ( l34[) and ( l35l) the equality of the put prices for the low strikes 

Vx < x , Vy < y* ai (x),P ai (x,y) = P a2 (x,y). 

For y > y^ix) = y& 2 (x), the equality is clear since both prices are equal to y — x. 

(vi) (i) : clear. □ 

Let us observe that the point (ii) of the last proposition allows to exhibit different 
volatility functions with analytic expressions that give the same put (or call) prices. Let 
us consider the same family as in subsection 14.21 coming from the call exercise boundary 
y{(x) = xf^. (assuming a,b,c > and max(c/a,6) < min(l, r/6)). For x > 0, we 
introduce the exercise boundary: 

yl{x) = y*(x) for x < x and y* 2 {x) = yl(x ) + (yl)'(x )(x - x ) for x > x . 

The condition (y% G C y ) is satisfied provided that (y*)'(x ) > max(l, £/r). This is au- 
tomatically ensured by the assumptions made on a, b, c since (y*)'(x ) = (bx^ + 2cx + 
ac)/{b 2 x1 + 2bcx + c 2 ). That family is such that cri(y) = d- 2 (y) for y < y^Xo). We can 

then calculate <t 2 as in subsection 14.21 using the relation aoix) = V 2 ^ ry 2^~ Sx ^y2( x ) j) rj^- 
gives (72 (x) = a(x) for x < x and for x > x , 



cr 2 (x) = 4/2 



(r(y*)'(x ) - 5)x + r(y$(x ) - x (yl)'(x ))] [((y$)'(x ) - l)x + y$(x ) - x (yt)'(x )] 



x 2 (yt)'(x ) 



In Figure [21 we have plotted the same example as in Figured] (x = 0.5 and y = 0.4), 
adding the graph of T t— > P a2 (T,x,y). The volatility function d2 has been calculated with 
the formula above with xq = 0.5. According to Proposition 15.21 and the Duality, the three 
prices are equal when T is large. In the second example (x = 3 and y = 1), we still observe 
that P a (T,x,y) and c^(T,y,x) converge toward the same value when T is large. On the 
contrary, the limit price of P U2 {T, x, y) is significantly different. To observe the same price, 
we should have taken, according to Proposition 15.21 x > 3. 



5.2 Calibration to "real" call and put prices 

In that subsection, we address some problems that arise if one tries to apply the calibration 
procedure when the prices p(K) and c(K) do not derive from a time-homogeneous model. 
We assume however that they are smooth functions of the strike K, and focus for example 
on the calibration to put prices. 
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Figure 2: P a (T,x,y) (solid line), P a2 (T, x,y) (dashed line with crosses) and c$(T,y,x) (dashed 
line) in function of the time T for a = 1, b = 0.4, c = 0.1, xo = 0.5, r = 0.2 and 5 = 0.1. 

Firstly, let us observe that the arbitrage-free theory allows to define a dual volatility as 
previously by (0, Y]: 



< Y, Vp (K) = -y . (36) 

Indeed, the payoff convexity in K ensures the positivity of p"{K) and the arbitrage-free 
assumption ensures that 5p(K) + K(r — 5)p(K) is nonnegative, so that the square- root is 
well defined. Let us prove the last point and suppose the contrary (i.e. 3y > such that 
^e Sy p{e^'^^ y ) < 0) to exhibit an arbitrage opportunity. In that case, there is z > y such 

that e Sy p(e {r ' 5)y ) > e Sz p(e ir ~ s)z ). We then sell one put with strike e ir ~ S)y and buy e 5{z ~ y) 
puts with strike e^ r ~^ z . This initial transaction generates a positive flow. The hedging 
works as follows: naming r the time at which the put sold is exercised, we have to pay 
e (r-s)y _ j n th. er words, we receive one share and borrow e^ r ~ s ^ y in cash. We keep this 
position until time r + z — y. At this time, we have exactly e s ( z ~ y ^ shares and puts with 
strike e^ r ~^ z . Thus, we obtain at least e s( - z ~ y ^e^ r ~ SS)z = e^ r ~ s ^ y e r ^ z ~ y ^ and we cancel the debt. 
The next proposition gives sufficient conditions that allow to construct an homogeneous 
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volatility which is consistent to the observed prices. 

Proposition 5.4. Let us assume that K e W + i— > p(K) is a C 1 function, C 2 on R"J_ — {Y} 
with Y = inf-fif > : p(K) = K — x } < +oo. Let us also assume that rj p defined by l[3b)) 
is bounded from below and above by two positive constants and admits a left-hand limit in 
Y. Then, if we extend r] p in any continuous function on (0, +oo) satisfying (H vo \) still 
denoted by r\ p , we have 

VK>0,P Vp {x o ,K)=p{K). 



Notice that once we choose the extended function r) p , we obtain rj p by first solving (1201) 
on starting from x with the condition y* p (xo) = Y and then using ( 1301) . 

Proof. The functions K t— > p(K) and K \— > c Vp (K,xq) solve (fT2"j) . Since we have < 
p{K) < K for arbitrage-free reasons, both functions go to when K — > 0. Thanks to 
Remark 11.21 they are proportional to #f and therefore there is A > such that: 

\/K<Y, p(K)=\c Vp (K,x ). 

The C 1 assumption made on p ensures p(Y) = Y — xq and p'iY) = 1. This gives 
gi(Y)/g'j(Y) = Y — xq and therefore Y = y* p (x ) using Lemma [3TT1 Thus, c Vp (Y, x ) = 
Y — Xq = p{Y) and A = 1. One concludes with Theorem 14.11 □ 

For the call case, everything works in the same manner, but we need to assume moreover 
that c(K) — > when K — > +oo. This is a rather natural hypothesis that plays the same 
role as p(K) — > when K — ► 0. 

Proposition 5.5. Le£ assume that K e i— > c(iT) zs a C 1 function, C 2 on — {X} 
TOt/j X = sup{K > 0, c(K) = Xq — X} > and limx^+oo — 0. i^ei us also assume 
that t] c defined by 



' /6V ; xy c"(x) 

bounded from below and above by two positive constants and admits a right-hand limit 
in X . Then, if we extend i] c in any continuous function on (0, +oo) satisfying (Tt vo i) still 
denoted by r\ c , we have 

VK>0,C Vc {x o ,K) = c{K) 
where rj c is obtained from rj c like a from a in the end of the proof of Proposition ^. 1\ 

Therefore, we are able to find volatility functions that give exactly the put prices and 
others that give exactly the call prices. Now, the natural question is whether one can find 
a volatility function a that is consistent to both the put and call prices. According to 
Proposition 15.21 all the volatility functions rj p (resp. r) c ) giving the put (resp. call) prices 
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coincide on (0, x ) (resp. (x , +00)). The only volatility function possibly giving both the 
put and call prices is 

{r] p (x) if x < xq 
7] C {X) it X > Xq 

We deduce from Proposition 15.21 : 

Proposition 5.6. Assume that ?] p (xq) = t] c (xq). Then, 

\/K > 0,p(K) = P a (x ,K) and c(K) = C a (x ,K) iff x*(Y) = x and T*(Jf) = x . 

6 The Black-Scholes model: the unique model invari- 
ant through this duality 

The purpose of that section is to put in evidence the particular role played by the Black- 
Scholes' model for the perpetual American call-put duality. We have recalled in the in- 
troduction that constant volatility functions are invariant by the duality. We have also 
mentioned that for the European case, the call-put duality holds for all maturities without 
any change of the volatility function. Here, on the contrary, we are going to prove that if 
the duality holds for the perpetual American options with the same volatility: 

Vac, y > P a (x, y) = c a (y, x) (37) 

then, under some technical assumptions, necessarily <t(.) is a constant function. 

Proposition 6.1. Let us consider a positive interest rate r and a nonnegative dividend 
rate 5 <r. We suppose that the volatility function a satisfies (7i V oi) ; o,nd is analytic in a 
neighborhood ofO, i.e. 

00 

3p > 0,Vx G [0,p), a{x) = ^a k x k . (38) 

Then, |3^ ) holds if and only ifWx > 0, cr(x) = a . 

We have already shown in the introduction that fl37|) holds in the Black-Scholes' case. 
So we only have to prove the necessary condition. We decompose the proof into the three 
following lemmas. 

Lemma 6.2. Let us consider a volatility function that satisfies (Ti, vo \). If the dual volatility 
function a is analytic in a neighborhood ofO, then the boundaries 2* andy*- are also analytic 
in a neighborhood ofO. 
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Lemma 6.3. Let us suppose that o satisfies (H vo i) and is analytic in a neighborhood of 
0. Let us assume moreover that r > 5. If the equality (Ulj) holds, a is constant in a 
neighborhood ofO: 

3p > 0,Vy G [0,p], a(y) = a . 
Lemma 6.4. Let us suppose that a is a constant function on [0, p] for p > satisfying 



(T~tvo\) a nd (37). Then, a is constant on M + (and x* and y% are linear functions). 



Proof of Lemma \6.°A Let us first show that x* is analytic in 0. Thanks to the relation 
we have p^J*^)) = y%{x) — x, and therefore = y — x*(y). Thus, x*(y) is analytic in 

iff (p{y) = is analytic in 0. Using the relation ( fl2l) and 0' = 1 — ^0, we get that is 
solution of 

<j>'{y) = l + J£—{{8- r)<j>{y)/y - 5{<j>{y)/yf) . (39) 

Notice that <ft(y) = y — x*(y) and ( 1T61) imply that if is analytic in then the coefficient 
of order in its expansion vanishes and the coefficient of order 1 belongs to (0, 1). 

To complete the proof we are first going to check that if ip(y) = YlkLi ( fiky k with 
0i G (0, 1) solves ( !39l) in a neighborhood of then = ip in this neighborhood. Then we 
will prove existence of such an analytic solution if). We have ip(0) = 0, and the function ijj 
being analytic with 0i ^ 0, its zeros are isolated points. There is therefore a neighborhood 
of 0, (0, 2e) where ip does not vanish. Let us consider 7 a solution of 7' — ^7 = starting 

from 7(e) 7^ in e : 7(0;) = 7(e) exp ^ -^jpjduj . Since ip solves (I3"9l . it is not hard to check 
that 7 is solution of ffT2l) with rj = a. The limit condition 7(0;) — * (cf. Remark II .2\ still 

x^O 

valid for ot when 5 = 0) is satisfied since we have -77-T ~ -r- and so f -rr^du — > —00. 

Thus we have 7(7/) = cg(y) with and ip(y) = g(y)/g'(y) = 4>{y) = y — x* a {y). We can 
then write x* a (y) = (f — 0i)w — ^^=2 * n ^ ne neighborhood of with 1 — 0i > 0. It is 
well-known that in that case, the reciprocal function y\ is also analytic in 0. 

Let us turn to the existence of ip. Since cr > q_ > 0, y — > is an analytic function 
in the neighborhood of 0. Thus, there is po > and a > such that 

2 00 00 
Vw G [0,p ], -r^p- = ^a k y k and ^ \ a k\Po < °°- 

° k=0 k=0 

The analytic function J2k>i ( t ) ky h solves (1391) if and only if 

00 00 / \ 00 / \ 

+ l)0fe+i/ = l + ((5-r)£ ai ^+ l - 6 J2[ Mi+i^J+i y k - 

k=0 k=0 \i+j=k / k=0 \i+j+l=k J 



Identifying the terms of order 0, we get that 0i solves -P(0i) = where P(x) = 5a x 2 + 
(I — (5 — r)ao)x — I. Since P{0) = — 1 < and P(l) = rao > 0, the polynomial P admits 
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a unique root on (0, 1) and we choose 0i equal to this root. Then, by identification of the 
terms with order k, we define the sequence ((pk)k>i inductively by 

^ i = (j - r ) E^=fc,^fc^0j+i - S E^+ i+ ;=fc,^fc,^fcQ^j+i^+i 

A; + 1 + (r — <5)a + 25a o 0i 

This ratio is well defined since (r — S)a + 2<5a o 0i = 5ao0i + l/0i — 1 > 0. 
We still have to check that the series Ylk>i 4>kV k is defined in a neighborhood of 0. To do 
so, we are going to show that there is p > such that the sequence (\4>k\p k )k>i is bounded. 
We have for 1 < k < n: 



k-l k ^ 

\S - r\J2\ak-j\p k ~ j \(f>j+i\p j + 

j=0 i=0 



E \0j+i\p j \4>i+i\p l | \ai\p* 

j-\-l=k—i 



Wi|p fc < 



k + 1 

Let us suppose that for 1 < j < k, |0j + i|p J < l/(j + 1). Then, 



k 



Wi|p K < 



k , *=0 \j+«=fc-i 



k + 1 

We remark that XI ]+ii+x = E jfi + T+i - 2 k=w^ > and we finall y S et: 

j+l=k—i j+l=k—i 

u fc< 2^ |M^ + p( |^- r | + 2 ^) Ej t 1 K-|^- 1 

l0fc+i|P < j^-q— j- (4UJ 

since ln ^~^ 1 2 ' >+1 < 1. Let us now consider A; such that VA; > fc , 2£|a | ln ^ 4 f 1 2 ' l+1 < 1/2. 
Now, we chose p G (0,p ) small enough such that VA; < k , \4>k+i\p k < l/{k + 1) and 
p(|<5 — r| + 25) yjjLi I^'Ip 7-1 < 1/2- Then we get by induction from (HOI) that VA; > 
fco,|0fc+i|p fc < l/(fc+l). □ 

Proof of Lemma 1 6'. ffi On the one hand, thanks to the assumption, a = a is analytic in 0, 
and therefore x* is analytic in thanks to Lemma [6.21 

oo oo 

3p > 0, Wy E [0,p),x* a (y) = ^xtf and a(y) = J^y*. 

i=l i=0 



On the other hand, it is not hard then to deduce from (1291) . a — a and the differential 
equation (|T9l that 

x:w= ^mp^Km . (4i) 
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From Corollary 13.61 and (I2"6"j) , we get 

(y ~ x*M)(ry - 5xl(y))((xl)- l )'(y) 



(M~ 1 (y)-y)(r(x*) 1 (y)-8y) 



(42) 



Now, we consider n = inf{« > 2,x« ^ 0} and suppose it finite. We can get easily that: 
x*(y) = x x y + x n y n + ... x* a (y)' = x x + nx n y n ~ x + ... 



and then 



(y) = £v- ~^y n + ■■■ iteWiv) = ±0- - if^ 1 ) + 



XI 



( 1 -K(y)/y)( r - Sx t(y)/y) 



(1 - x x )(r - 6xx) + x n (28x 1 - (r + 8))y n ~ L + 



x 



(y) 



s 



.1 i 



X T 



y j \ v 

The right hand side of (1421) has then the following expansion: 



(1 - x x )(r - 8x1) + — ((r + 8)x x - 2r)y 



n-l 



,-y ' ■ 



+ . . . 



H 1 + 



(1 — x x )(r — 8x1) 



25xi - r + 8 + — - -^j- 

X n 



n-l 



nx 



+ 



The equality of the terms of order n — 1 in fl42|) then leads to: 

x n 



n—l 

TlX yi.X 1 



(1 — x\){r — 8x\) 
Since x n 7^ and with a simplification we get 

1 



[28x^ +1 - (r + - (r + + 2r] 



nx r 



nil + x^ 1 ) 



r — 8x1 



n-l 



28x1 + (r - 8)J2 x i+ 2r 



k=l 



(43) 



In the case 5 = this gives n(l + x™ x ) = x™ 1 + • • • + xi + 2 which is not possible because 
Xi G (0, 1). When < 8 < r, we denote a = r/8 > 1 and rewrite (|43l : 



n(l + Xi )(ct-xi) = -2x™ + (a-l)x"~ i H h(a-l)xi + 2a = ct-x™ + (a-xi 



1 - x? 
1 — Xi 



, l-z? 



Therefore, n(l + x? 1 ) = Xl + ^— — < 2^—^- because /3 1— > 4-^- is decreasing on (1, a) 



n—l fe 
X . 



(x™ < Xi). To show that this is impossible, we consider P n (x) = n(l + x n 1 ) — 2 ^2 k=0 
WehaveP„,(l) = and for x < 1, P' n (x) = n(n- l)x n ~ 2 -2 J2lZl kx k ~ l = 2 A;(x"- 2 - 



x 



< 0. Thus P n is positive on [0, 1) and P n {x\) > which is a contradiction. 



Proof of Lemma \6.4\ It is easy to get from (TT^j) and = that 

x *Mv{<{y)) 



<(y)' 



y°{y) 



□ 



(44) 
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We have a(x) = cr for x G [0,p]. Since x*(y) solves (jHJ) and x*(y) < y, x* a {y)' = x*(y)/y 
on [0,p\. Therefore, x*(y) = X\y for y G [0,p\. Thanks to (fTUj) . X\ is the unique root in 
(0, min(l, r/6)) of 

xictq = 2(1 — a;i)(r — 5x\). 

Now let us observe that (HU) gives for y G (0,y|(p)], ac*(y)' = 2 (y-^(y) ( Kry-fa»( y )) with 
= Xip. Since y — » a?iy solves this ODE, for which local uniqueness holds thanks 
to the Cauchy Lipschitz theorem, we then have x*(y) = X\y on [p, y%{p)] and so y%{p) = 
(x*) _1 (p) = p/xi. Then, gives a /a(y) = 1 on [p, p/xi]. Thus, we prove by induction 
on n that x* (?/) = a?iy and cr(y) = o"o for y G [0, p/ {x\) n }. This shows the desired result. □ 

7 Conclusions and further developments 

Addressing Call-Put duality for American options with finite maturity in models with time- 
dependent local volatility functions like ([T]) would be of great interest. For the perpetual 
case treated in this paper, we could take advantage of a very nice feature : in the continu- 
ation region, the price of the option writes as the product of a function of the underlying 
spot price by another function of the strike price. Unfortunately, this product property no 
longer holds in the general case. 

Next, according to our numerical experiments (see figure [2]), American Put and Call 
prices computed in infinite maturity dual models may differ for finite maturities. This 
means that in the case of a time-homogeneous primal local volatility function q(t, x) = cr(x), 
if there exists a dual local volatility function for some finite maturity T, then this volatility 
function is either time-dependent or depends on the maturity T. On the contrary, in the 
European case presented in the introduction, time-homogeneous volatility functions are 
preserved by the duality. 

Let us nevertheless conclude on an encouraging remark. Let P(T, x, y) denote the 
initial price of the American Put option with maturity T and strike y in the model ([T]) and 
x*(T,y) stand for the corresponding exercise boundary such that P(T,x,y) = (y — x) + if 
and only if x < x*(t,y). Then the smooth- fit principle writes 

P{T,x*{T,y),y) = y-x*{T,y) 
d x P(T,x*(T,y),y) = -l 

Differentiating the former equality with respect to y yields 

d x P(T, x*(T, y),y)d y x*(T, y) + d y P(T, x*(T, y),y) = l-d y x*(T,y). 

With the second equality, one deduces that d y P(T, x*(T, y),y) = 1. Therefore the smooth- 
fit principle automatically holds for the dual Call option if there exists any. 
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